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CN 1 Abstract 



CO • We propose a set of nonlinear integral equations to describe on the excited states of an 

integrable the spin 1 chain with anisotropy. The scaling dimensions, evaluated numerically 
in previous studies, are recovered analytically by using the equations. This result may be 
relevant to the study on the supersymmetric sine-Gordon model. 
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CN 
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tJ- ; 1 Introduction 
O 

The ID spin systems have been providing problems of both physical and mathematical interests. 
See e.g., pJEl- Among them, there exists a family of solvable models of the Heisenberg's type 
with spin-S 1 [31 Q]- In this report, we are interested in the excited states of a member in the 
c~| , family, the spin 1 chain with anisotropic interaction. 

The recent progress in the study of the integrable system brings forth a powerful machinery, 
the method of the nonlinear integral equations (NLIE, for short) j^EI- The NLIE method has 
been successfully applied to the study of the XXZ model (S = ^). It clarifies the finite size 
property of the ground state as well as the scaling behavior of excited states (of corresponding 
six vertex model (Zj). The application is not restricted to the lattice models: it also provides the 
detailed descriptions of the excited states in the field theoretical models such as the sine-Gordon 
model El EDI CUE] and perturbed conformal field theories [EH IbTj. 

The study on the higher spin cases has, however, encountered technical difficulties. This has 
been resolved in [lj, at least for the ground state. There the NLIE, which is relevant to the 
evaluation of the free energy for arbitrary S, is derived for the isotropic case . See also ^HI f° r 
an interesting application of the result to the 0(4) nonlinear sigma model in the limit S — > oo . 

In this report, we extend the study to excited states of the anisotropic S = 1 chain . Simple 
assumptions, suggested by numerical investigations, lead to a set of NLIE which enables the 
evaluation of energies for arbitrary system size. The proposed NLIE has a structure which seems 
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to be a natural extension of the excited NLIE for the spin 77 case. We will analytically verify the 
previous observations on some low lying excitations by numerical methods |171 118j . 

The result obtained here may be not only relevant in the spin chain problem. Recently the 
study on the excited states in supersymmetric sine-Gordon model attracts much attentions |19l 
I2()| . It is expected that the proper discretization of the model is given by the inhomogeneous 
anisotropic spin 1 chain. We thus hope that the current study will shed some light on the analysis 
of the supersymmetric sine-Gordon model. 



2 The model and the assumptions 

We are interested in the spin 1 chain with anisotropic interactions. The hamiltonian contains 
several two-body interactions, 



H=J2(^ - (^) 2 + cos2 7 K z - (af) 2 ) - (2 cos 7 - + ~ 4sin 2 7 (Sf) 2 ) (1) 

i=i 

where a short-hand notation 

Ui = Si - S i+ i =a i + of 

is employed. For simplicity, we impose the periodic boundary condition, S^ +1 = Sf, a = x,y, z. 
The parameter 7 specifies the anisotropy. Throughout this report, we only consider the range 
7 < § and exclude 7 of the form 7 = 2£ with (m, n) co-prime. 

We follow the strategy in |15j and start from a integrable 19 vertex model [211 122 j . The 19 
vertex model is obtained from the 6 vertex model by the fusion procedure [22j 122]. The latter 
one is associated to the spin ^ XXZ model while the former corresponds to the spin 1 chain. 
The strategy is to treat the spin 1 model and the spin ^ model simultaneously. To be more 
precise, we introduce commuting transfer matrices T\{x) and T2(x) acting on spin 1 quantum 
space consisting of ./V sites. The auxiliary space for the former one is given by the spin \ space 
while it is spin 1 for the latter. 

Their explicit eigenvalues read, 



Ti(z) = <p\z - i7j — — — \-<p(z + t"f) 



T 2 {z) = cP(z - i-)(p(z - 1—)—- — + <P(z - i-)(p(z + z-) 



Q{z) " v " Q{z) 
.37 g^+jg) 

2'^ % 2 ' Q(z-q) +nZ "2*^ ' u 2' Q(z + q)Q(z-q) 

2'^ 1 " 2 ' Q( z + q) 



sinh^ z 



Note that in place of standard spectral parameter u , we choose u = ix . The important function 
Q is given by the Bethe ansatz roots Zj, (j = 1, • • ■ , M); Q(z) = YlfLi sinh(z — Zj). We will 
denote three terms of T%(z) in © by Aj(x), i = 1, 2, 3. 

The eigenvalue E of the Hamiltonian is evaluated though 



L ^IAi og r 2 (x)U_ i7/2 . (3) 
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Thus once if Q is obtained, the evaluation of E is straightforward. This is equivalent to say 
that, if all locations of the Bethe ansatz roots are known for excited states, then the energy is 
evaluated. This is, however, a formidable task for large systems. Apart from few lower excitations, 
it is extremely difficult to find all locations of Bethe ansatz roots. 

The most crucial observation in the NLIE formulation is that this task is avoidable. With 
proper choice of auxiliary functions, one can bypass dealing with a complete set of Bethe ansatz 
roots; one only has to deal with finite number of complex roots characterizing the excitation 24 , 
171 11UI ITS1 114j . This may break down if very high excitations are of our interest. We nevertheless 
believe that this method will be efficient for the treatment on excited states which are physically 
important in the thermodynamic limit. 

Let us be more accurate. In the ground state, the Bethe ansatz roots are given by the sea of 
2 strings. By a 2 string we mean a pair of Bethe ansatz roots Xi ± iyj with \yj — ^| -C 1. 
We consider excitations for which only a finite number of roots deviate from the sea of 2 strings. 

In addition, we will make the following three assumptions. 

Assumption 1 

There should not be a pair of complex roots z\, Z2 such that Z\ — Z2 is a multiple of 7. 

For the string type solutions, it is known that the separation of neighboring roots in a string 
deviates slightly from 7 for finite system sizes [2S]. Our assumption thus does not contradict 
with this pattern. It, however, excludes the complete strings 2B1 127j. Therefore we should devise 
another route to deal with the case when q(= e* 7 ) is at roots of unity, which is beyond the 
present scope. 

Assumption 2 

The following classification of complex roots , other than 2 strings, are possible. 

1. inner roots : \$Xzj\ < Z, j = 1, ••■ , Mj 

2. close roots: ^ < \Qzj\ < j = 1, • • • , Mq 

3. wide roots: ^ < \^szj\ < = 1, ■■■ ,Mw 

4. self conjugate roots : \3>Zj\ = f , j = 1, • • • ,M p . 

The above classification has already been proposed in which discusses the excitations in 
the limit TV" — > 00. There a complex conjugate pair of inner roots is referred to as a narrow pair 
while a pair of close roots is to a intermediate. We adopt a different notation for similar roles 
played by them in comparison with the spin ^ case. We will sometime denote the locations of 
inner roots by s - , close roots by , wide pairs byu^ 1 and self-conjugate roots by pj + ^i. Here 
the upper index +(— ) means its imaginary part being positive (negative). 

The zeros of T\{x) and T2{x) play also important parts. The numerical investigation for N 
up to 8 leads to the remarkable feature. 

Assumption 3 

The zeros of T\(x) and T%(x) in the strip Qx S [—7/2,7/2] distribute exactly on the real axis. 
We denote their locations by 9j, (J = 1, • • ■ , N\) and by 9a \ (« = !,-•• , iY 2 ) , respectively. 



3 



For example, we plot the zeros of T\{x) and T2(x) for two cases in fig.^ The state in the left 
figure in ^ corresponds to an excited state in the singlet sector (M = 6), system size is 1 N = 6 
chain and with the coupling constant 7 = 4^. For the state in the right figure inQ parameters 
are chosen M = 7, N = 8 and 7 = . The unit of tics in imaginary direction is normalized 
to \. These figures thus clearly support the assumption3, which simplifies the derivation of the 
nonlinear integral equations drastically. 



Figure 1: Zeros of T\{x) and Tiix), Left: an excited state in iV = 6 chain M = 6, 7 = ^% Right: 
N = 8 chain M = 7(S Z = 1), 7 = ^. 

We supplement the explicit locations of corresponding BAE roots. 



BAE roots for JV = 6, M = 6, 7 = ^ 


0.65640208387 +1.5707963268 i 
-0.11533405190+0.20504539087 i 
-0.24721889635 


-0.53549516834 +1.5707963268i 
-0.11533405190-0.20504539087i 
0.35698008462 



BAE roots for N = 8, M = 7, 7 = ^ 


0.61918637956+ 1.5707963268i 
-0.83691028264+ 1.5707963268i 

-0.18745636322 
0.29966102387+ 0.15214900339 i 


-0.10487718381+ 1.5707963268i 
0.29966102387 -0.152149003388i 
0.024363369664 



3 Auxiliary functions and Sum rules 

In this section, we introduce several auxiliary functions which are crucial for our purpose. 
Firstly we define the most natural auxiliary function a(z), defined by 

_ \ 2 {z + i%) _ X 3 (z - if) 

a[z) ~ \i(z + q) ~ x 2 (z-q) 

In view of a(z), the Bethe ansatz equation can be cast into the form a(zj) = — 1, (j = 
1, • • • ,M). One also uses its logarithmic form, Z^\zj) = 2irlj , where a(z) = exp(iZjp (z)), 
which leads to the root density function formulation in the thermodynamic limit. The first 
auxiliary function, a(z), thus has the deep connection to the Bethe ansatz equation and plays 
the fundamental role in the study of the spin ^ chain. The previous study shows that, 
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unexpectedly, this is not the case for general values of S at their ground states. Instead, the 
most crucial functions for the ground state for 5 = 1 are given by 



:= Ai(x) + A 2 (x) - := A 2 (x) + A 3 (x) 
X s (x) Ai(ar) 

Physically, b\(x) is related to the density function associated to the centers of 2-strings. 
For a technical reason, we introduce the shifted functions, 1 

b(x) = b\(x + ie) b(x) = b\{x — ie) and capital ones 25(x) = 1 + b(x), 2$(z) = 1 + b(x). The 
definition obviously concludes, 

T 2 (x) = 4>{x + il)4>{x + - ie) (5) 



2' rv 2'Q^+n) 

2 ' a'o^-f) 



^-4) 9 , ( r + Y? ^+^) (6) 



(2) 

and 53 (x) = when x = - — ie. 

In analogy to a{z) and Z^^x) we introduce zff(x) := Alogb(x). In contrast with 

zff (x),zff(x) is in general a complex- valued function. We assume that its real part is an 
almost monotonic increasing function of x and that the imaginary part vanishes when the real 
part takes integers (half- integers). One then applies the similar argument for the spin ^ ^01 to 
derive the following sum rule, 

N 2 = 2(S + M p ) - -#(3 7 S + 0(2 7 S)) + 2M W + M c (7) 

where 9{x) := tt[^ + ^\ and [^J specifies the integer part of x. 

The validity of this rule is checked against many examples. The rule is crucial in the deter- 
mination of the constant term in the NLIE. 

We have a remark. The repeated integers are observed for zff(x), which are attributed to 
the existence of special holes/roots ^3 Our case studies indicate no symptom of "special 
holes/roots" for Z N (x). We therefore dismiss the possibility in this report. Even if they exist, 
only a small modification will be needed in the following argument. 

We need to introduce another pair of auxiliary functions, one of which coincides with T 2 , 
apart from the normalization. 

yi{x):= Ji -37W Y^x) := 1 + yi (x) (8) 

<p(x — i-f)4>{x + l-f) 

The essential observation to derive the NLIE for In b remains almost same as the one made 
in the case of the largest eigenvalue sector of the quantum transfer matrix [Td]. 
We introduce renormalized functions, 

T 2 V (x) = r 2 (x) := ]J tanh f (x - 9^) (9) 

r 2 (x) - LJ - 27 



1 In case of the finite temperature problem, further shifts in x are needed for analyticity reason, which is not 
necessary for the finite size problem. 
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and consider the integral, 



-ikz 



dz, 



where C encircles the strip of width 7 (FigEJ) in the counterclockwise manner 





y/2i 


1 
1 

-t 




1 
1 
[ 


i 




-y/2i 



Figure 2: The integration contour C. The straight line in the lower half plane is termed as Cg 
while the one in the upper plane is referred to as C u after reversing the direction. 

Thanks to the renormalization factor, T% (x) contains no zeros or poles inside C. The Cauchy's 
theorem thus concludes 







-ikz 



dx 



c, 



-ikz 



dx. 



(10) 



where Cg, C u means the lower and the upper part of the contour, respectively. 

One then substitutes T% in terms of 25 for C u and in terms of 25 for Cg by utilizing 
) and (0). This provides the various relations among auxiliary functions in the Fourier space. 
After straightforward manipulations we find the desired NLIE 



In b(x) 



C b + iD b (x + ie) 
+ I G s {x-x')\n^{x')dx' 

7 



+ 



-00 
00 



K(x - x' 



G s (x - x' + 2ie) ]n<B(x')dx' 
i + ie)lnY 1 (x')dx' 



11) 



In yx{x) 



D y (x) + 



[l\(,r - ./ + |i - ie) ln*8(z') + K(x - x' - -U + ie) \D.*&(x'))dx' 

(12) 
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where kernel functions read 
G s (x) 



K(x) = — 



1 r°° sinh(f - ^)k 
2^ 2cosh(^A;)sinh(f - 7 )fc 

1_ f°° 1 
2~ 



e ikx dk 



Jkx dk. 



, 2cosh(^A:) 

It is worth mentioning that G s (x) is related to the logarithmic derivative of the soliton- 
soliton scattering matrix of the sine-Gordan model, and K(x) is a standard kernel function in 
the thermodynamic Bethe ansatz equation for the RSOS model [2^ 
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The constant C& is found by matching both sides of NLIE at x — > — oo and it reads, 

C b = mS -iOinS). (13) 

We remark that, depending on choice of branches, this value is determined only modulo 27ri. 
This ambiguity can be absorbed into the definitions of branch cut integers. 

The driving term in (fTT|) consists of three parts: D&(ar) = ^buik( x ) + ^hcL( x ) + ^roots( x )- The 
first term is identical to the one for the ground state, thus we refer it to as the bulk contribution, 

^buikC*) = N Xk{x) 

where Xk( x ) = 27rK(x) and xk(— °o) = 0. The second consists of two pieces, the contributions 
of holes of T\ and T2 . 

bSLw = !>*(*- of) + E^- ^ 2) ) 

j a 

where x( x ) is an °dd primitive of 2ttG s (x) and x(0) = 0. The third term represents the contri- 
butions from complex zeros other than 2 strings, 

dSL(x) = - 53 X ( x -(c°-aii))+ 53 U ~ {w {~ a¥ \ + Y.^^r-)^) 



v=±,j <r=±,3 



... 1, sinh(x — id) 
4>s(x) := - log 



(15) 



also be written as 



i sinh(x + id) 

77(77—37) 

2(tt-27) 



and parameters are given by 5 = if^zr^ an d ^ = 1 — IT- The last two summations in (j!5[) can 



where we adopt a notation; for any function f(x), fn(x) = f(x) + f(x — i~/sgn(^sx)) ^ 

The driving term for lnyi lacks the bulk contribution and composed of two terms, D y (x) = 
^ol + ^roL- Explicitly, 

D^M^^XKix-e^ + li) Di y Ux) = -i J] aXK ( x - c ° +ai i). 

a o-=±,j 

For given locations of excited zeros and holes, eqs. (|11|) and ()12j) fix the values of auxiliary 
functions (modulo exp(2ivi)). Then the evaluation of energy spectra for any iV is immediate by 
the following expression, 

£ = e + e ho ie + e roo ts + ee. (16) 

In the above eo denotes the bulk ground state energy eo = —N(cotj + cot 27). The second 
term ehoie stands for the excitation energy for a hole, 



7T 



7 cosh -x 

' 7 
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Among the contributions from complex excitations, the one from the close roots appears here 
explicitly, 

e r oots = - / J e(cj - ia-). 

<7=±,j 

This phenomenon is also observed for the spin ^ case ^U] for 7 < ^. 

The last term contains implicit contributions from all excitations, and it is given by 

&b = — J K\x — x' — ie + i^)\ x _,_ il i2^^{x')dx' + - j K'{x — x'+ie— i^-) |a; >~ i-y/2 hi 2$ (a/) da;'. 

The locations of complex roots and holes are, however, not given a priori. Therefore an improve- 
ment of the set of NLIE is necessary so as to make the evaluation of these locations possible. To 
resolve this, we need the NLIE for In a with a general complex argument z. 

The derivation of the NLIE can be done in two ways. One starts from the derivation for 
real z, then apply the analytic continuation argument in JU] to derive the equation valid for y 
with larger complex parts. The procedure is straightforward but for a small technical difficultly 
which does not show up in the spin ^ problem. Alternatively, one can start directly from y with 
larger complex part, apply carefully the following simple lemma ^ to reach desired NLIE. 

Lemma 1 We define, for a smooth function f(z), its "Fourier transformation" by 

1 f'°° 

fvik} = ^ J f{x + iy)e- lkx dx. (17) 

with z = x + iy. 

Suppose f(z) has a pole at z = zq with the residue r and analytic elsewhere. For Qzq > 
then 

' f[k]e~ ky 0<y<Qz 
(f[k] - ire~ ikzo )e- ky y > %z 



fy[k] 

Similarly for Qzo < 

' f[k]e~ ky $z <y<0 



fy[k] 



(f[k] + ire~ ikz °)e~ ky y<%z 



We have derived the following equations by these two manners and checked that they lead to 
identical results. 

A moment consideration convinces us that we need to treat the equations, at least by three 
separate regimes for positive imaginary values of z. 
For the simplest case, Qz S [0, ^) , the NLIE reads 

/OO /'OO 
K{z-x -ie)\n^dx'+ J K(z - x' + ie)\n¥>dx (18) 
-00 J —00 

where the constant and drive term are given by 

C a = 2i9(2jS) - i6(3jS + 6(2jS)) 
D a {z) = J2 XK(z-(c°-ali))-J2xK(z-e&). 

a=±,j a 



s 



The equation ceases to be valid as crosses ^ due to the singularities existing in the kernel 
functions of the above equation. Taking account of the modification due to them, we find the 
equation for G , ^) 

lna(z) " ln ( 1 + ^3-)) = -^-^-^) 

- ! G s {z-i^ -x' -ie)\n<2>dx' + ! G s {z - £ - x' + it) In <&dx' 

./-oo ^ ./-oo ^ 

/oo 
^(z-x')ln(l + yi(x'))^' (19) 
-oo 

The simultaneous evaluation of a(z) at different strips is thus necessary within this framework. 
When z = ry, some spurious singularities appear in the rhs, which causes some numerical 
problems. We hope to report on the remedy in a separate communication. 
The last case, Qz <G § ), the equation takes again a simpler form, 



G^iz-ij- -x' -ie)lnmx' + / (2 - - x' + ie) In %>dx' (20) 

-oo ^ ./-oo ^ 

Note that integration constant is null modulo 2iri and the drive term reads, 

D a (z) = 2^^( )- 2^ M „ ) 

a ' c=±,j 

_ E tf( ^-q-K-.f.) ) _ Etf(i -q- j , ) 

cr=±,j ^ ^ 

Although these three NLIE for In a take rather involved forms, their meaning is simple: once 
In 53, In *B and In Y are known on the real axis, the NLIE yield the evaluation of In a at arbitrary 
z. Again, we have checked the validity of these NLIE at various points in the complex plane. 

As noted above, the locations of complex roots and holes are determined by "quantization" 
conditions 

In b(^ 2 )-ie) = (2/^ + 1)™ 
]na(wf) = (2/| w±) + l)vri, 

(21) 

To be precise the NLIE still leaves 2-rri ambiguity as remarked in the above . Thus one has 
to be careful in the choice of branch cuts. We will present several examples of proper choices in 
the next section. 

The set of properly chosen integers then fixes the locations of complex roots and holes, 
thereby the energy spectra. In this sense, the NLIE characterizes the finite size spectra com- 
pletely. 



lna(0f j ) = (2/j aj + l)7d, 
lna(c0 = (2/] c±) + l)m, 
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4 Low lying Excitations in the thermodynamic limit 



For an illustration, we consider the energy levels of some low lying excitations in the thermody- 
namic limit. We consider the weak anisotropy 7 — > case. As shown in I24| ITU] . the scaling 
behavior of these levels can be evaluated without solving the NLIE explicitly. 

The contribution to the excited spectra mainly comes from the left and right extremes of 
roots distribution. In the limit, N — > 00, these two contributions are almost decoupled. 

The energy scales as 

E ~ Ne + -^A£ 

where eo denotes the bulk ground state energy. The scaling energy AE is given by the sum of 
left and right contributions, AE = AE + + A EL. We conveniently introduce scaling functions, 

a±(x) = a(±-(x + IniV)), A±{x) = 1 + a±(x) 

7T 

b±(x) = b(±-(x + lniV)), B±{x) = l + b±(x) 

7T 

y±{x) = yi(±-(x + IniV)), Y±(x) = 1 + y±(x) 

IT 

Similarly, to holes near the left (right) extremes, we associate the new locations 9 + (9~), 

ef = ±^ef-io g N (22) 

Then the left/right contributions are explicitly written as 

AE± = —Y ± — / e- {x±ier) log B±dx T — [ e~^ e ') log B±dx (23) 
7 ■ *7 J l l J 

where e' = — e. In the following discussion, it is irrelevant thus we set e' — > for notational 
simplicity. 

The NLIE also can be transformed into scaling forms. We prepare associated kernel functions 
and drive functions, 

K~'{x) = -K(^x), G"i{x) = -G s (^x) 

7T 7T 7T 7T 

X l {x) := x{-x), Xk(x) ■= Xk(-x) 

Then the set of resultant NLIE takes the form, 

lnb±(x) = ^2ie~ x + C± + iD±{x) + J CT{x - x') ]nB±(x')dx' - J G^(x - x') ]nB±(x')dx' 

+ J K^{x - x' ^^i)lnY±(x')dx' 
lny±(x) =C± + D±(x) + J K^(x - x' ± i|) hxB ± {x')dx' 

+ J K 7 (x-x'=F^)ln B ± {x')dx'. 
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Those scaling NLIE for lna± have similar but more involved forms. We only write the case 

lna±(x) = C± + iD± - J K^(x - x') log B±(x')dx' + J K^(x - x') log B ± {x')dx' 

The constants and D functions in the above depend on the choice of specific state of interest 
and choice of branch cut integers. Below we consider simple examples for which the results of 
numerical studies are available. 
Example 1: 

We consider the lowest excitation in the sector S z = 1. It is shown numerically that the N — 2 
zeros of Q(x) form the sea of the two strings and the last zero is located at the origin[T7|ll8j. In 
the strip Qx € [— ^, T\ possesses no zeros while two zeros of T2 are located on the real axis. 
We denote their scaled locations (|22|) by 0^. 
The drive terms in this case are found to be, 

C± = T(*-S)h D± = ±x 1 (x-6 ± ) 

C± = ±ni, £>± = ±ix y K (x -9 ± ± 

We conveniently choose the quantization conditions, log b±{9 ± ) = =p(2/ ± + l)-7r. The locations 
of holes then satisfy 



2i 



^-(271^ + 5)i- J G^(9 ± -x) In B±(x)dx 
+ J G^(9 ± - x) In B ± (x)dx - J K" l {0 ± - x T |i) In y±(x)dx) (24) 

Our choice of and of the branch cut integers leads to I 1 * 1 > in the present case. By 
substituting (|24j) into (|23|). we present AE± only in terms of auxiliary functions, 



AE± = ^-(± J e~ x log B±dx± J e~ x \ogB±dx + (27r 2 / ± + 5n)i 

± J ' K^{9 ± -x±^i)\nY±{x)dx} (25) 
The standard dilogarithm trick leads to the explicit values of desired integrals, 
4i(± j e~ x \og B±dx =F J e~ x log B±dx 

^Jg^-^b^^Jg^-^b^ 

±tt J K' y (9 ± -x=p ^») In Y±(x)dxj 
± (tt — 2<5)i(log B±(oo) — logi?±(oo)) =F i7rlog Y±(cxd) 
= j \\ogb±)'\ogB±dx - j \ogb±(\ogB±)'dx + j (log b±)' log £±cte - J log 6±(log B±)'dx 

+ J (log J/±)' log Y ± dx - J (log y ± =F 7rz)(log Y±)'dx (26) 
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where we have used x^(oo) = tt It is worth mentioning the asymptotic values 

6 ± (-oo)=0 6±(-oo) = 

y±(-oo) = -1 

1 _|_ e T^h _ i _|_ e ± 2 «7 

b±(oo) = 6±(oo) = e=F2 , 7 

y±(oo) = 1 + 2 cos 27. 

Using these values and by the change of integration variables, we find that the rhs of (|26|) is 
given by the dilogarithm functions, 

2L+(6±(oo)) + 2L+(6 ± (oo)) + 2L + (y±(oo)) - 2L+(1) + 2L(1) + 2L + (1) T 7rilogF±(oo) (27) 

Explicit definitions are as follows. 

t i \ 1 r^ 1 °g( 1 + y) l °gy ^ 

L ^ = 2h { — TTy 1 

m := -i /*(!S«Pzrt + **V )dy 
2 Jo y 1 - y 

The scaling energy AE is divided into two parts AE = AE\ + AE2. The former brings the 
central charge AE\ = while the second is related to the scaling dimension, AE2 = 2irvX. 
The sound velocity is readily evaluated v = — . 

It is easily shown that the sum of the first four terms in (|27jl remains constant for small 7. 
Thus the dilogarithm formula, utilized in the study of the ground state JH] in the rational case 
(7 — > 0), is also applicable. We identify these terms as the contribution to AEi = ^c. This lead 
to the correct central charge |. 

Evaluating the remaining contributions, we find, 

2ir 2 1 
AE 2 = (/ + + /_ + X c + -) 

7 o 

where X c = '■ ■ Therefore, by choosing the minimal value /+ = 1_ = 0, the above calculation 
recovers the desired scaling dimension, X = X c + g [T71 ITS] . 
Example 2: 

Let us consider another excitation in the sector S z = 1. The N — 2 zeros of Q function again 
form the sea of the two strings while the last zero is located at ^i. In this case, T\(x) possesses 2 
zeros on the real axis while T 2 (x) does 4 zeros. We denote corresponding scaling locations by ^ 
for zeros of T\(x) and Of ,{j = 1> 2) for those of T 2 (x). In this case it is convenient to introduce 

Xk( x ) := Xk{x) - it so that Xk(°°) = °- 
Then the drive terms are then given by, 

3 



Ct = D} = ±i ~ Of ± |i) 



3 



ct = £>r = =F2^x^-^ 

J 
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The quantization conditions read 



logb ± (9f) = T (2/f + 1)tt loga ± (e ± ) = ±(2.7* + l)vr 
One easily verifies Ij , J ± > 0. 

The quantization condition for logfr 1 * 1 leads to an expression analogous to (|24|) . This time, 
however, the rhs contains ^ Y2j Xjc^j ~ The quantization condition for logo 1 * 1 enables us 
to represent this by integrals, 



J2^ 6 f-^ = - 7T± - [ K^(t-x)logB ± (x)dxT- [ Ki{t-x)\ogB±{x)dx (28) 
. i J i J 

We used a simple relation x~k( x ) + xtc(~ x ) = ~~ 71 an< ^ P u * ^ = ®- 
We then proceed as example 1 and obtain the same AE± and 

2tt 2 , 7 v 2vr 2 . 1 1, 

7 7r(7r — 27J 7 7 loA c 2' 

where we choose I-=j — l. This again coincides with the numerical result |17j . 

We also analyzed several examples in the different spin sectors, and checked that the results 
all recovered the desired values. We, however omit further discussions for brevity. 



5 Discussion and Summary 

In this report, we derive a set of NLIE which characterizes excited state spectra of the spin 1 
chain with anisotropic interactions. The equations are tested numerically against exact data. 
Some desired scaling dimensions are derived analytically for some low excitations. 

Finally, we comment on an implication of the result obtained here. Through the light cone 
approach |30| 131 j . the inhomogeneous version of the spin 1 chain, or the 19 vertex model will be 
the proper candidate for the discrete analogue of the supersymmetric sine-Gordon model . The 
latter's excited spectra have been attracted attentions recently |20| . A proper deformation of 
the NLIE obtained in this paper may be useful in such investigations. Assume that the several 
conjectures on the analytic properties in the spin chain problem are also valid for supersymmetric 
sine-Gordon model. Then it is readily shown that we reach the almost same nonlinear integral 
equations. The "bulk" driving term of log b should be then replaced by mL sinh ^x where mL is 
related to the inhomogeneity A by mL = 2A"exp(— ^A). We however avoid drawing a conclusion 
in haste as it requires careful analytic and numerical checks; the consistency to the S— matrix 
picture in |22j, for instance. We hope to report this in a near future, together with complete 
discussion on the conformal limit. 
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